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Abstract 

Time-Delay Interferometry (TDI) is the data processing technique needed for generating inter- 

o 

cr*j ' ferometric combinations of data measured by the multiple Doppler readouts available onboard the 

three LISA spacecraft. Within the space of all possible interferometric combinations TDI can 
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generate, we have derived a specific combination that has zero-response to the gravitational wave 
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fact it works for waveforms of any kind. Our approach is analogous to the data analysis method in- 



signal, and called it the Zero-Signal Solution (ZSS). This is a two-parameter family of linear com- 
binations of the generators of the TDI space, and its response to a gravitational wave becomes null 
when these two parameters coincide with the values of the angles of the source location in the sky. 
Remarkably, the ZSS does not rely on any assumptions about the gravitational waveform, and in 



troduced by Giirsel & Tinto in the context of networks of Earth-based, wide-band, interferometric 



gravitational wave detectors observing in coincidence a gravitational wave burst. The ZSS should 
be regarded as an application of the Giirsel & Tinto method to the LISA data. 
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I. INTRODUCTION 



LISA, the Laser Interferometer Space Antenna, is a deep-space mission of three spacecraft 
flying in a triangular formation. Jointly proposed to the National Aeronautics and Space 
Administration (NASA) and the European Space Agency (ESA), it aims to detect and study 
gravitational radiation in the millihertz frequency band via laser interferometry [lj]. 

Modeling each spacecraft as carrying lasers, beam splitters, photo-detectors^and drag- 
ree proof masses on each of two optical benches, it has been shown jjj, 0], p, Q, 
Bl> Q that the six measured time series of Doppler shifts of the one-way laser beams 
between spacecraft pairs, and the six measured shifts between adjacent optical benches 
on each spacecraft, can be combined, with suitable time delays, to cancel the otherwise 
overwhelming frequency fluctuations of the lasers (Aw/v ~ 10~ 13 /VHz), and the noise due 
to the mechanical vibrations of the optical benches (which could be as large as Avjv ~ 
10 _16 /\/Hz). The achievable strain sensitivity level h ~ 10 _21 /v / Hz is set by the buffeting 
of the drag-free proof masses inside each optical bench, and by the shot noise at the photo- 
detectors. 

In contrast to Earth-based, equal-arm interferometers for gravitational wave detection, 
LISA will have multiple readouts. The data they generate, when properly time shifted 
and linearly combined, provide observables that are not only insensitive to laser frequency 
fluctuations and optical bench motions, but at the same time show different couplings to 
gravitational radiation and to the remaining system noises. The technique for synthesizing 
all these interferometric combinations has been called Time-Delay Interferometry (TDI), 
and it has been shown that the functional space it generates can be obtained by prop erly 
time-shifting and linearly combining the four Sagnac generators (a, 0, 7, £) 4(, jl^ . 
In particular, there exists a combination of the four generators that, for a given waveform 
and source location in the sky, achieves maximum signal-to- noise ratio (SNR) Q|. 

Here we further elaborate on the geometrical properties of the four generators' responses. 
We find that there exists a two-parameter family of combinations of (a,/3, 7, C) that has 
null response to an observed gravitational wave signal. In the case of gravitational wave 
bursts, or sinusoidal signals observed over a time scale short enough that the amplitude and 
frequency modulation induced on the received signal by the LISA motion around the Sun are 
negligible, this new TDI combination shows a sharp null when its two parameters coincide 
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with the values of the angles describing the source location in the sky, or its mirror image 
with respect to the LISA plane. This combination, which we call the Zero-Signal Solution 
(ZSS), does not rely on any assumptions about the gravitational waveform, and in fact it 
works for gravitational waveforms of any kind. Although its ability to identify these two 
points in the sky degrades at lower frequencies, it should be regarded as an improvement 
over the Symmetric Sagnac combination, £, |l2| since it could be used for discriminating 
gravitational wave signals from spurious fluctuations generated by instrumental noise in the 
entire LISA frequency band. 

The ability of the ZSS to identify the source location in the sky follows from the fact that 
the time spent by the gravitational wavefront to propagate across the LISA constellation 
enters into the gravitational wave responses of the four generators (a,/3, 7, £)• The two 
independent travel times spent by the wave to propagate from one spacecraft to the other 
two define the "six-pulse" U structure in the responses of the four generators. Since the 
signal contribution to the ZSS combination becomes null at the correct values of these two 
time-delays, and the antenna pattern of the ZSS is symmetric with respect to the LISA 
plane, it follows that in general there are two possible points in the sky where the signal 
must have originated. If the duration of the signal is such that the effects of the motion 
of the LISA array around the Sun can be disregarded (as in the case of bursts), then the 
two-fold degeneracy cannot be removed. However, if the duration of the signal is such that 
the effects of the motion of LISA around the solar system barycenter can be measured (as in 
the case of signals from binary systems), then the source location can be uniquely identified. 
As LISA moves around the Sun, the point in the sky where the source is located will remain 
fixed, while the other will move in such a way to remain its mirror image with respect to 
the LISA plane. By constructing the ZSS with different stretches of data, the unique source 
location can be identified. 

The ability of the ZSS to identify two possible points in the sky where an observed 
gravitational wave signal must have come from will be particularly useful when observing 
gravitational wave bursts, or continuous signals that are strong enough to be observed by 
integrating the data over time scales shorter than the time scale during which the phase 
and amplitude modulations of the signal become measurable. In the latter case, the ZSS 
should be regarded as the first step of a hierarchical procedure for measuring the parameters 
of a continuous signal. The ZSS could be particularly useful when observing, for instance, 



3 



gravitational radiation emitted during the inspiral of two super-massive blackholes, or one 
smaller object into a super-massive hole. The ability of the ZSS to identify the location 
of the binary system in the sky, without relying on the knowledge of the signal waveform, 
might reduce significantly the number of templates needed in order to study waveforms from 
such astrophysical systems [111 ]. 

This paper is organized as follows. Section II provides a brief overview of Time-Delay 
Interferometry, the data analysis technique needed to remove the frequency fluctuations of 
the six lasers and other noise sources affecting the LISA measurements. After showing that 
the entire space of interferometric combinations can be generated by properly combining the 
"six-pulse" Sagnac combinations (a, /3,j, £), in Section III we obtain the general expression 
of the ZSS in the Fourier domain, and note that it does not rely on any assumptions about the 
time-dependence of the observed gravitational wave signal. The corresponding expression 
of the ZSS in the time-domain is then given in Section IV, where we show that the three 
data sets, a, /3, and 7, enter with eighty-one different time delays. Section V derives the 
expression of the ZSS in the long-wavelength approximation (i.e. when the characteristic 
wavelength of the wave is much larger than the typical linear size of the LISA array). In 
Section VI we give a summary of our work and emphasize that this "part I" article only 
addresses the geometrical properties of the ZSS. A quantitative analysis covering its ability 
to discriminate gravitational wave signals from random noise fluctuations affecting the LISA 
data, and the accuracies it will achieve in reconstructing the source location in the sky and 
the wave's two independent amplitudes, will be presented in a companion article. 



II. SUMMARY OF TIME-DELAY INTERFEROMETRY FOR LISA 



A brief overview of TDI, the technique needed to remove the frequency fluctuations of the 
lasers from phase measurements obtained with an unequal-arm interferometer, is presented 
in this section. The reader is referred to the artie.es Q.l . Q. B, Q for a derivation of TDI 



valid for the non-rotating and rigid 
(TDI1)), and to the following ones 



SA configuration (the so called "first generation TDI" 



2 1, P Q for the modified expressions accounting 
for the LISA spacecraft motions (the "second generation" TDI (TDI2)). The overall LISA 
geometry is shown in Figure 1. There are six beams exchanged between the LISA spacecraft, 
with the six Doppler measurements j/y = 1, 2, 3) recorded when each transmitted beam 
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FIG. 1: Schematic LISA configuration. Each spacecraft is equidistant from point o, with unit 
vectors pi indicating directions to the three spacecraft. Unit vectors hi point between spacecraft 
pairs with the indicated orientation. 

is mixed with the laser light of the receiving optical bench. The frequency fluctuations from 
the six lasers, which enter in each of the six Doppler measurements, need to be suppressed 
;o a level smaller than that identified by the secondary (proof mass and optical path) noises 
l| in order to detect and study gravitational radiation at the predicted amplitudes. 

Since the LISA triangular array has systematic motions, the two one-way light times 
between any spacecraft pair are not the same . Delay times for travel between the space- 
craft must be accounted for depending on the sense of light propagation along each link 
when combining these data as a consequence of the rotation of the array. Following [10], 
the arms are labeled with single numbers given by the opposite spacecraft; e.g., arm 2 (or 
2') is opposite spacecraft 2, where primed delays are used to distinguish light-times taken 
in the counter-clockwise sense and unprimed delays for the clockwise light times (see Figure 
2). Also the following labeling convention of the Doppler data will be used. Explicitly: 
JJ23 is the one-way Doppler shift measured at spacecraft 3, coming from spacecraft 2, along 
arm 1. Similarly, y 32 is the Doppler shift measured on arrival at spacecraft 2 along arm 1' 
of a signal transmitted from spacecraft 3. Due to the relative motion between spacecraft, 
L\ 7^ L\ in general. As in pj], we denote six further data streams, Zij = 1,2,3), as 
the intra-spacecraft metrology data used to monitor the motion of the two optical benches 
and the relative phase fluctuations of the two lasers on each of the three spacecraft. Finally, 
the light times and L i (i = 1,2,3) are not only different (pure rotation) but also change 
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FIG. 2: Schematic diagram of LISA configurations involving six laser beams. Optical path delays 
taken in the counter-clockwise sense are denoted with a prime, while unprimed delays are in the 
clockwise sense. 

in time. For this reason, in the subscript notation for delays, their order is important jsj] , Q] , 
[lfl| . The frequency fluctuations introduced by the gravitational wave signal, the lasers, the 
optical benches, the proof masses, the fiber optics, and the measurement itself at the photo- 
detector (shot noise) enter into the Doppler observables jjij, with specific time signatures. 
They have been derived in the literature ^], Q], and we refer the reader to those papers for a 
detailed discussion. The Doppler data y^, are the fundamental measurements needed to 
synthesize all the interferometric observables unaffected by laser and optical bench noises. 

Let us consider the TDI2 Sagnac observables: (pj Ji P2,P3) 3|- Their expressions in terms 
of the Doppler measurements yij, Zij are as follows 



Pi = [V31 + 2/23:2 + |/12;12 + 2/21;312 + 2/32;3'312 + J/l3;l'3'312] 

— [2/21 + V32;3' + 2/l3;l'3' + 2/31;2'l'3' + 2/23;22'l'3' + 2/l2;122'l'3' 
+ -[(^21 — £3l) — (Z21 — ^3l);2'l'3'312 + (^32 — ^12);3' 

+ (^32 — ^12j;12 — (^32 — ^12);3'312 — (^32 — ^12) ;122'1'3' 
+ (^13 ~~ z 23);2 + (^13 — ^23j;l'3' ~~ \ Z 13 — ^23);22'1'3' 

— (^23 — ^13);1'3'312] , 



(1) 



with p 2 , P3 following from Eq. (JTJ by permutations of the spacecraft indices. The semicolon 
notation shown in equation ((TJ emphasizes that the operation of sequentially applying two 
or more delays to a given measurement is non-commutative, and a specific order has to be 



adopted to adequately suppress the laser noises 3, jl(3]. Specifically: jjij-ki = Hij(t — Li(t) — 
L k (t - Li)) ^ i/ij-ik. 

The expressions of the gravitational wave signal and the secondary noise sources entering 
into pi will in general be different from those entering into a, the corresponding Sagnac 
observable derived under the assumption of a stationary LISA array 

a, a. 

However, 

the other remaining, secondary noises in LISA are so much smaller, and the rotation and 
systematic velocities in LISA are so intrinsically small, that index permutation may still 
be done for them llOf. It is therefore easy to derive the following relationship between the 

□ n 

signal and secondary noises in px, and those entering into the TDI1 combination a p, plfl 

P x{t) ~ a(t)-a(t -L x - L 2 -L 3 ) , (2) 

where L, , i = 1,2,3 are the unequal- arm lengths of the stationary LISA array Equation 
(J2J) implies that any data analysis procedure and algorithm that will be implemented for 
the TDK combinations can actually be derived by considering the corresponding TDI1 
combinations. For this reason, from now on we will focus our attention on the gravitational 
wave responses of the TDI! obse^ables ( a , A7 , C ) S, fl. 

In the reference frame co-moving with LISA and centered at the point o that is equidistant 
from the three spacecraft (see Figure 3), the response of a to a gravitational wave signal is 
given by the following expression 3 



gw 



1 - -^(/i 2 - /is) (*i(t - I - L x - L 2 ) - tf x (t - /x 3 I - L 2 )) 



+ 



1 + — (// 2 - JU 3 ) 

1 - -p-Oa -Hi) 

Li 2 
I / 

1 + —{Hs-Hi) 

J->2 



+ 



/ 

T 3 
i 



(/ii - /i 2 ) 



1 + -r(Hl- H2) 



($i(t - // 3 I - Lx - L 3 ) - - ii 2 I - L 3 )) 

(* 2 (t - us i - l 2 ) - * 2 (t - m /)) 

(^ 2 (t - m I - L 2 - L x - L 3 ) - * 2 (i -Hsl-Lx- L 3 )) • 
(tf 3 (t - // x Z - L 3 - La - L 2 ) - * 3 (t -ii 2 l-Lx- L 2 )) 
(V 3 (t-ti 2 l-L 3 )-* 3 (t-iixl)) ■ (3) 



Here pi = k ■ pi, with j3j , i = 1, 2, 3 being the unit vector pointing from o to spacecraft i, 
and k the unit propagation vector of the wave. The function is equal to 

1 hi ■ h(t) ■ hi 



2 i-ik-hiY 



(4) 
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FIG. 3: The array and the wave coordinates. The x, y, and z axes are the LISA coordinate axes 
with the origin at the point o, and the array lies in the x, y plane. The x', y', z' axes are the 
coordinate axes for the incoming gravitational wave. The z' axis is parallel to the direction of 
propagation of the wave. The x', y' axes are those with respect to which h+(t) and h x (t) are 
defined. The angles 9, <ft describe the location of the source in the sky. 



where h(t) is the first order spatial metric perturbation at point o. Note that L\k ■ fi\ = 
1(1^2 — ^3), and so forth by cyclic permutation of the indices. The gravitational wave h(t) can 
be written as [h+(t) e + + h x (t) e x ], where the 3-tensors e + and e x are transverse to k and 
traceless. With respect to an orthonormal propagation frame (x',y',z') their components 
are equal to 



( 1 N 
0-10 

\0 J 



f 1 N 

1 
v o j 



(5) 



Equation (jSJ) for a gw shows a "six-pulse response" to gravitational radiation. In other words, 
a 5-function gravitational wave signal produces six distinct pulses in (a, /3, 7) which are 
located with relative times depending on the arrival direction of the wave and the detector 
configuration. 

Together with (a, j3, 7), the symmetric Sagnac combination, (, has been shown to form a 
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basis for the TDI1 space of combinations and, like a, has a six- pulse response to gravitational 
radiation. Furthermore, it can be written in terms of (a, (3, 7) in the following way jjj 

C - C.123 = « 1 - <*, 23 + Al _ A23 + 7,1 - 7,23 • (6) 

The four interferometric combinations (ct,/3^j,C) jointly also give the expressions of the 
interferometric combinations derived in j^: the Unequal-arm Michelson (X, Y, Z), the 
Beacon (P, Q,R), the Monitor (E, F, G), and the Relay (U, V, W) responses 



Xi 


= a 23 _ A2 - 7,3 + C , 


(7) 


p 


= C - a 1 , 


(8) 


E 


= a - C,i , 


(9) 


U 


= 1,1-13. 


(10) 



The remaining expressions can of course be obtained from equations (0 |Hl El nTJ|) by permu- 
tation of the spacecraft indices. All these interferometric combinations have been shown to 
add robustness to the mission with respect to failures of subsystems, and potential design, 
implementation, or cost advantages 

a, a. 

One important property of the £ combination is that it is much less sensitive to 
gravitational waves than any of the other Sagnac combinations (a, [3, 7) in the low part of 
the frequency band accessible to LISA (/ < 5 mHz). Since, however, it is affected by the 
same instrumental noise sources, it was recognized that it can be used as a "gravitational 
wave shield" for estimating the spectra of the instrumental noises affecting the "sensitive" 
combinations l^j], 3]. 

Although the ( combination will be a very useful tool for assessing the on-fhght-noise 
performance of LISA , , the question of whether there exists a TDI combination that 
completely suppresses the gravitational wave signal still remains open. In what follows we 
will tackle and solve this problem. 

III. THE ZSS IN THE FOURIER DOMAIN 

LISA has the capability to simultaneously observe a gravitational wave signal with several 
interferometric combinations of different antenna patterns and noise transfer functions j^J. 
For this reason it should be regarded as an overlapping array of interferometric gravitational 



wave detectors working in coincidence. For a given signal waveform and source location in 
the sky, it has in fact been shown Q| that there exists an optimal combination of the 
generators of the TDI space that achieves maximum SNR. 

I I 

In what follows we will answer a complementary question to the one addressed in jl4j . 
We will identify, within the LISA TDI functional space, a specific interferometric combina- 
tion that has zero-response (and therefore zero SNR) to a gravitational wave signal. Such 
a combination provides a way for discriminating a gravitational wave signal from noise- 
generated fluctuations, and for this reason it should be regarded as an improvement over 
the Symmetric Sagnac combination £ jl^ . 

Our approach is closely related to the one introduced by Giirsel and Tinto |l3| in the con- 
text of data analysis for networks of ground-based interferometers observing in coincidence 
a gravitational wave burst. There, a network of three ground-based, wide-band, detectors 
widely separated on the Earth was shown to provide enough information for uniquely solving 
the "Inverse Problem" in gravitational wave astronomy: the reconstruction of the wave's 
two independent amplitudes and the determination of the location of the wave source in 
the sky. With three non-coplanar detector responses and two independent time delays, it 
is possible to identify a unique, two-parameter linear combination of the three independent 
detector measurements that has zero-response to the gravitational wave signal when the two 
par ameters coincided with the values of the polar angles of the source location in the sky 

For LISA there exists an analogous two-parameter linear combination of the four gen- 
erators of the TDI space, (a,/3, 7, £), which in the high-frequency part of the LISA band 
is capable of identifying two points in the sky where the observed gravitational wave burst 
must have come from. The property of this combination, which we have named the Zero- 
Signal Solution, follows from the fact that the two independent travel times spent by the 
burst propagating across the three spacecraft define the separation of the "six-pulses" in the 

n 

responses of the four generators (a, /3, 7, Q |2j|. Each travel time defines a circle on the sky 
where the signal must have come from, and their intersection reduces the source location to 
two possible points. However, due to the co-planarity of the antenna patterns of the four 
generators, the LISA ZSS is unable to resolve the two-solution ambiguity for bursts, as was 
possible with networks of three Earth-based laser interferometers 131 ] . 

In order to derive the expression of the ZSS, let us consider in the Fourier domain an 
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arbitrary element of the TDI space, which we will refer to as r\ 

fj(f) = ai(f, A) a(f) + a 2 (f,X) P(f) + a 3 (f, A) j(f) + a 4 (/,A) £(/) , (11) 

where the {a,i(f, A)}| =1 are arbitrary complex functions of the Fourier frequency /, and of a 
vector A containing parameters characterizing the gravitational wave signal (source location 
in the sky, waveform parameters, etc.) and the noises affecting the four responses (noise 
levels, their correlations, etc.). For a given choice of the four functions {a>i}f = i, rj gives an 
element of the functional space of interferometric combinations generated by (at, (3, 7, £). 
The corresponding expression of the signal-to-noise ratio, SNR"^, of the combination rj is 
equal to 

SNRl = f fu |Ql 5 * + a ^ + «3 % + aZ\ 2 df (12) 
Jfi (\ai a n + a 2 f3 n + a 3 7 n + a 4 Cn| 2 ) 

In equation (JT2J) the subscript n refers to the noise part of (5, (5, 7, Q, the angle brackets 

represent noise ensemble averages, and the interval of integration (fi, f u ) corresponds to the 

frequency band accessible by LISA. As a consequence of the relationship between ( and a, 

(3, and 7 (equation ©), it has been shown (see equation (10) of reference |lj]) that equation 

|T2|) can be rewritten in the following equivalent form 

snb> = i°.g-+*fr-+«.yi' df . (13) 

Jfi (|ai a n + a 2 p n + a 3 7„| 2 ) 

For a given choice of the three functions {aj}f =1 , SNR^ gives the corresponding SNR of an 
element of the TDI space, whose Fourier transform is equal to rj = a\ a + a 2 (3 + 03 7. Our 
goal is therefore to identify, for a given gravitational wave signal, the three functions {aj}f =1 
that give SNR^ = 0. Since SNR^ is a positive definite functional, it can be zero if and only 
if the numerator of its integrand is identically null, i.e. 

fj^(f) = a x aP°{f) + a 2 &»{f) + a 3 i 9W (f) = . (14) 

The expressions for a gw , (3 9W and 'y 9W can be rewritten in the following form (see equations 
(J3J), (HJ), and appendix B) 

a 9W (f) = a + (f,£>s,cj> s ) h+(f) + a x (fAM hx(f) , (15) 

p gw (f) = MM.) h + (f) + MMs) h x (f) , (is) 
7H/) = -y+ifAM K(f) + 7x hM) , (17) 
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where (0 S , S ) are the two angles describing the location of the source in the sky with respect 
to a LISA coordinate frame (see Figure 3), and the functions a+, a x , /3 + , /3 X , 7+, 7x are 
pattern functions, derived in appendix B. From these considerations it follows that rf w can 
be rewritten in the following form 

fT(f) = K a + (f,e s ,<f) s ) + a* MAM + a 3 MAM) M) 

+ [ai a x (f,9„<p.) + a 2 P x (f,6 s ,(f> s ) + a 3 7x (/, 9 S , S )] K (f) ■ (18) 

which is obtained by substituting equations (|T51 ITol 117)1 into equation (f*HJ) . The wave's two 
independent amplitudes, h + , h x , are referred to the wave axes (x',y') which, without loss 
of generality, can be assumed to be oriented in such a way that the x'-axis lies parallel to 
the LISA plane (the x — y plane). In order to have rj 9W = for any arbitrary pairs of wave's 
amplitudes h + , h x , equation (fTHj) implies that the following homogeneous linear system of 
two equations in three unknowns must be satisfied 

ai a + {f,e s ,(f> s ) + a 2 (3 + (f,8s,<Ps) + a 3 7+ (/,0„0.) = (19) 
ai a x (f,9 s ,(f) s ) + a 2 (3 x {f,9 s ,(f) s ) + a 3 7x (/, 9 8 , (j) s ) = (20) 

Since the rank of the matrix associated with the linear system above is in general equal to 
two, it is easy to derive the expressions for a\, a 2 , a 3 that give zero- response to a gravitational 
wave signal at the source location 

ai (f,9 s ,(f> s ) = p + (f,0sAs) ixifAM - MAM i+UAM > (21) 

a2(f,9 s ,<f> a ) = l+(f,9 s ,(f) s ) a x (f,9 s ,(f) s ) - 7x (/,0 a ,&) a + (f,9 s ,(j) s ) , (22) 

a 3 (f,e.M = MAM M,0sAs) - MAM MAM ■ ( 23 ) 

Equations (|211 12*2*1 12^|) imply that the following linear combination of the Fourier transforms 
of the three generators (a,/3,j) 

fj= \p+tfA<t>)MA<i>) - MA4)MA<i>)]c*(f) 
+ [ 1+ UA<f>)MA<l>) - MAtiMAWHf) 
+ [ a+ (/, 0,0) MA<!>) - MA 4) MA 4)] lU) , (24) 

is a function of the two parameters (9, 0), and it has null response to the gravitational wave 
signal observed in (a,/3,j) when (0,0) — > (0 S ,0 S ), regardless of the particular waveform 
considered. 
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Note that the expression of rf w given in equation ()24|) changes sign under reflection with 
respect to the LISA plane. This is because the functions a+, 7+ are invariant under 
reflection with respect to the LISA plane, while a x , (3 X , j x change sign under the same 
operation (see appendices A and B for the derivation of these symmetry properties). This 
implies that the combinations a + (3 x — a x /3 + , /3+7 x — /3 X J+, 7+a x — 7x« + , which define the 
ZSS fj, also change sign under reflection with respect to the LISA plane. In the case of a 
gravitational wave burst, or if we assume to observe a periodic signal for time scales shorter 
than the time during which the Doppler effect induced by the LISA motion around the Sun 
begins to be observable, we conclude that the ZSS will not be able to distinguish between 
the source location in the sky and the point that is its mirror image with respect to the 
LISA plane. 

In the specific case of a sinusoidal signal, for which the two independent polarization 
components are characterized in the source rest reference frame by a constant frequency 
fo, it is interesting to calculate the time scale during which such a periodic signal will still 
appear in the LISA data as a sinusoid. In the LISA reference frame, any of the TDI responses 
will measure a sinusoidal signal whose frequency is Doppler modulated by the LISA motion 
around the Sun in the following way 

/(*) = /o[l - k -^\ , (25) 
c 



where k is the unit vector of propagation of the wave, V is the velocity of the LISA guiding 
center o relative to the solar system barycenter, and c is the speed of light jl7l |. If we 
take the root-mean-square (r.m.s.) value of the frequency change given by equation (J25j) . 
averaged over all possible incoming directions of the GW signal, we obtain the following 
resulting r.m.s. value of the Doppler change of the signal frequency in the LISA reference 
frame 

A/=((/(t)-/ ) 2 ) 1 / 2 = ±/o^, (26) 

where Vq is the velocity amplitude of the almost sinusoidal motion of the array around the 
Sun, and the angle brackets () denote an ensemble average over source directions. If we 
take Vq ~ 30 km/s, and assume fo = 10~ 2 Hz, the corresponding r.m.s. frequency shift is 
equal to Af = 5.0 x 10~ 7 Hz. In other words, if we would integrate the data for about 20 
days, the data frequency resolution would be equal to or larger than the corresponding r.m.s. 
frequency shift induced by the Doppler modulation. Since the frequency shift induced by the 
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amplitude modulation is smaller than the Doppler frequency modulation at this frequency 
fo, ( 3] '13)' it follows that over a period of about 20 days, the expression of the ZSS can 
be used for identifying two points in the sky where an observed sinusoidal signal might have 
arrived from. 

As an example, in Figure 4a we plot contours of constant power spectral density of the 
ZSS in terms of its two-parameters, (9,4>), in the case of a circularly polarized sinusoidal 
signal of frequency fo = 10~ 2 Hz, incoming from 9 S = 32°, S = 107° and for noiseless 
data. Note that the ZSS is equal to zero only at the correct source location, and at its 
mirror image with respect to the LISA plane. Also, the plane of symmetry (defined by 
6 = 90°) shows an identically null ZSS. This is because every TDI response is insensitive 
to the "x" polarization components of signals incoming from sources located in the LISA 
plane (see Appendix A). This makes the ZSS degenerate for these source directions. Figure 
4b shows instead the power spectral density of the ZSS as a function of 9, for the cases 
= (f) s = 107°, and = 207°, while in Figure 4c we plot it as a function of and for the 
choices 9 = 9 S = 32°, and 9 = 59°. From these two plots it is possible to have an idea of 
the sharpness and narrowness of the nulls in terms of the angles (9, <fi). In order to estimate 
the magnitude of the accuracy by which these two points in the sky can be identified with 
this method, however, it will be necessary to include the effects induced by the noises and 
the magnitude of the SNR by which the signal is detected in the three Sagnac observables 
(a,/3, 7) [13]. This analysis will be presented in our forthcoming "part II" article. 



IV. THE ZSS IN THE TIME DOMAIN 



The expression of the ZSS can also be given in the time-domain by inverse Fourier trans- 
forming equation (j24j) . Having the ZSS in the time domain will be particularly useful when 
attempting to discriminate gravitational wave bursts against statistically significant fluctu- 
ations generated by instrumental noise affecting the three TDI combinations (a, (3, 7). A 
quantitative analysis of the improvement in statistical confidence added by using the ZSS 
as a veto against noise fluctuations will be addressed in a follow up paper. 

In order to derive the time-domain expression of the ZSS, let us consider only one of 
the three terms entering into the expression of the Fourier transform of the ZSS given in 
equation (J2I|), say as(f, 9,<p) 7(/). This is because the other two terms entering into the 
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FIG. 4: (a) Contours of constant level of the modulus-squared of the ZSS constructed for a sinu- 
soidal signal of frequency fo = 10~ 2 Hz, incoming from a source located at: 9 S = 32° , 4> s = 107°. 
The plot is symmetric with respect to the LISA plane (9 = 90°); (b) two cross sections along 
the #-axis of the modulus-squared of the ZSS for the same signal of (a); one curve corresponds 
to 4> = <p s = 107°, and the other to <j> = 207° / 4> s (c) two cross sections along the 0-axis, with 
6 = e s = 32°, and 6 = 59°. 
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ZSS can be obtained from it by permuting the spacecraft indices. 

Since the function as(f, 9, <fi) can be written in the following analytic form 

a3(/,^0) = E4 fc) (^0)e 2 ^ A 3 fc) , (27) 
k=i 

where the twenty seven "amplitudes" and the corresponding "delays" A^ depend only 
on the angles (9, <p) and not on the Fourier frequency /, it follows from the linearity of the 
inverse Fourier transform that the expression of the ZSS in the time domain assumes the 
following form 

27 

V (t, 9, 0) = J2l4 k) <*(t - A[ k) ) + 4 k) P(t - ) + 4 k) 7 (f - A?>)] . (28) 
k=i 

If we take into account that each of the Sagnac observable shows a gravitational wave pulse 
at six-distinct times, we conclude that for an arbitrary pair (9,(f>), the ZSS will display 
in principle "four hundred eighty six" pulses. At the correct source location, and at its 
mirror image with respect to the LISA plane, the ZSS becomes null, making all these pulses 
disappear. 

Figures (5, 6, 7) schematically show the implementation of the ZSS in the case of a 
gravitational wave burst whose wave components, h + (t), h x (t), are delta-like pulses of unit 
amplitudes. The noiseless responses of the three Sagnac combinations to such a gravitational 
wave signal, which we have assumed to be incoming from (9 S = 32°, <p s = 107°), are given 
in Figure 5 (a, b, c). Note the characteristic six-pulse structure of their responses. In 
Figures 6 (a, b, c) the time-domain expressions of the inverse Fourier transforms of the 
functions ai(/, 9, <p) &(/), a 2 (/, 9, <f>) /3(f), a 3 (f, 9, 4>) j(f) are plotted for values of the angles 
{9 = 31.5°, = 106.5°). The corresponding ZSS combination is given in Figure 6 (d), 
showing that no cancellation is achieved for values of the polar angles slightly different from 
those corresponding to the source location. Figure 7 (a, b, c, d) shows instead plots similar 
to those given by Figure 6, but now with (9, 4>) = (9 S , <j> s ). 

V. THE LONG- WAVELENGTH LIMIT 

In order to understand the behavior of the ZSS, 77, in the low-part of the LISA frequency 
band, we provide in this section its analytic expression in the long- wavelength limit (LWL), 
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FIG. 5: The responses of the three Sagnac combinations, a (a), (3 (b), and 7(c), to a gravitational 
wave burst of unit amplitudes. The source has been assumed to be located at: 9 S = 32° , 4> s = 107°. 
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FIG. 6: The time domain responses of the Inve^e-Fourier Transforms of the three functions, a± 
(a), a<2 (3 (b), as 7 (c), which determine the ZSS. The angles (9,4>), over which the search 
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FIG. 7: Plots (a, b, c, d) are as in Figure 6, but now with (0 = 6 S , 



i.e. when the characteristic wavelength of the GW signal, A, is much larger than the typical 
scale size of the array, L. 

Let us rewrite the expression for rj in the Fourier domain (see equation (|24j) above) 

rj = [0+(/.M7x(/,M - /W, M t+(/. Ml 
+ [t+C/.MW/.M) - 7x(f,o,<P)<x + (fA<t>)]P(f) 

+ M/,M) ftc(/>M) - «x(/,M) /?+(/> M)] t(/) , (29) 

where the functions a+, a x , fix, and 7+, 7 X are given in appendix B. We have been able 
to find the long-wavelength expansion (f'L << 1) of rj by using the program Mathematica 
[lsI ]. Since the LISA arm lengths will differ at most by a few percents the general 
expression of the ZSS in the long wavelength limit can be written in the following simpler 
form 



fj(f, 9, 0) « (2 tt/ I/) 4 [(F+F 2 X - F*F+) + (F+F 3 x - F 2 X F 3 +) + (F 3 +F* - F s * J? 

x 



4 

+ri 3 ■ hW (/) ■ n 3 ] + a n + fi n + %} , (30) 

where L is the nominal LISA arm length (~ 5 x 10 6 km), the label ^ on a specific functions 
means "evaluated at 6 = 9 S ; <fi = <p s ", and (a n , (3 n ,j n ) are the Fourier transforms of the 
three random processes associated with the noises in (a,/?, 7). Note that as (k x , k y , ±k z ) — > 
(k x ,k y ,±k z Y s \ rj — > as expected. This is because the unit vector k, corresponding to a 
point that is the mirror image of the source location relative to the LISA plane, will have 
the same scalar product with the three unit vectors hj as k^. 

It is interesting to compare the LWL expression of the ZSS given above, against that for 
the Symmetric Sagnac combination, ( |l2| 

C 9W (f) * \~ 2 (2^/L) 3 [(* ■ AOCna • h(/) • h x ) + (k ■ n 2 )(h 2 ■ £(/) ■ h 2 ) 

+ (k - n 3 )(n 3 • h(/) • n 3 ) + - a n + /?n + 7n , (31) 

where the noise affecting £ can be rewritten (in the LWL) in terms of the (a, /3, 7) noises 
by using the LWL limit of the Fourier transform of equation (jHJ) with equal arm lengths. 
Although both combinations display the (27r/L) 3 -dependence in the signal and the same 
noises, the ability of rj to further suppress the gravitational wave signal by identifying the 
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direction it came from should be regarded as an improvement over (. This property will 
be particularly useful for discriminating exceptionally strong signals against random noise 
fluctuations that will be still observable by ( in this low-part of the LISA frequency band. 

VI. SUMMARY AND CONCLUSIONS 

We have identified a method for nulling a gravitational wave signal observed in coincidence 
by three LISA TDI combinations. We have named the resulting TDI combination the 
Zero-Signal Solution. The method can in principle identify two points in the sky where a 
gravitational wave burst, or a sinusoidal signal, came from. In order for the method to work 
in the case of a sinusoidal signal, the gravitational wave signal must be observable after 
integrating the data over a period such that the Doppler modulation induced by the LISA 
motion around the Sun can be disregarded. 

Although the analytic expression of the ZSS has been derived in terms of the three Sagnac 
observables (a, /3, 7), any other three TDI combinations capable of generating the entire TDI 
space (such as the Unequal- Arm Michelson combinations X, Y, Z) could have been used. In 
particular, the TDI combinations (A,E,T) Q], whose noises are uncorrelated, would have 
worked as well. The use of the (A, E, T) combinations will be particularly useful when we 
will estimate the accuracies provided by the ZSS in determining the source location and the 
wave's two amplitudes. These points will be addressed in our follow up article. 
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APPENDIX A: THE BEAM PATTERN FUNCTIONS 

This appendix derives the analytic expressions of the antenna patterns Fj~(9, <j>), F* (9, (j>) 
entering in the wave functions ^j(t) , j = 1,2,3, and in the pattern functions a+, a x , P+, 
j3 x , 7+, 7x introduced in equations (fT51 ITBl IT7|) . and given in appendix B. 
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Consider a plane gravitational wave of direction of propagation k incident on the LISA 
spacecraft constellation. Let the wave coordinate system be (x',y',z'), with the wave trav- 
eling in the z' direction, and let the axes with respect to which h + (t), and h x (t) are referred 
to be the (x',y') axes. The coordinate system associated with LISA, (x,y,z) is obtained 
from them by a rotation described by the polar angles (9, <fi) (Figure 3). For sake of sim- 
plicity, and without lack of generality, we can assume the x' axis to be parallel to the (x, y) 
plane, and the wave's two independent amplitudes h + (t), and h x (t), associated with the two 
independent polarizations, to be defined with respect to this choice of coordinates. 

Let us now consider the null vector rh defined by 

rh = -^{zx' + i <v) , (Al) 

where e x >, e y r are the unit vectors in the x', y' directions, respectively. The tensor h is then 
just H 

h = 2 h+(t) Re{rh <g> rh) + 2 h x (t) Im(rh <g> rh) . (A2) 

In the LISA coordinate system the components of hij can be obtained by finding the com- 
ponents of rh in these coordinates. The vector fh can be written as 



1 

V2 



rh = — -j= [(— sin (f> + i cos 9 cos (f>) e x + (cos (ft + % cos 9 sin 0) e y + {i sin 9) e z ) , (A3) 



while the unit vector k, describing the direction of propagation of the wave from the source 
to LISA assumes the following form 

k — (— sin 9 cos <p) e x + (— sin 9 sin <p) e y + (— cos 9) e z . (A4) 

Since the ^j(t) , j = 1,2,3 functions, entering into the expressions of a 9W , /3 gw , j 9 ™ have 
the following analytic forms 

T . . 1 hi • h(t) - hi . . . 

#,(t) = - -1 — — I , A5 

after replacing equation (jA2|) into (|A5|1 . and performing some simple algebra, we can rewrite 
*&j{t) in the following form 

* (t) s F n d ^) h + (t) + F*(9,<p) K(t) 
2{l-{k.h j f) 
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where the functions Fj~(9, 0), F*(9, 0) are equal to 

Ff(6, 0) = 2 [fij ■ Re{m <g> rh) ■ fij) , (A7) 
Ff(0, 4>) = 2 [nj • Jm(m <g> m) • , (A8) 

If we now rewrite the vector rh in the following form 

fh = fhi + i rfi2 , (A9) 

where the two vectors rh\, rfi2 have components 

rhx = -j= [(-sin0) + (cos0) e y ) , 

m 2 = — ^= [(cos 9 cos0) + (cos# sin0) e y + (sin#) e z ] . (A10) 
v2 

the tensor product rh®fh becomes equal to 

rh ® m = [(mi (g) To]) — (m 2 (g ^2)] + * [(^1 <8> ^2) + (^2 <8> ^1)] • (AH) 

The expressions for Fj~, F* given in equations (|A7[ IA8jl can be rewritten as 

Ff = 2 [(fij ■ rh,) 2 - (fij ■ m 2 ) 2 ] , (A12) 
F/ = 4 (% • mi) (hj ■ m 2 ) , j = 1, 2, 3 . (A13) 

Since the unit vector hj are in the (x, y)-plane, they can be expressed as 

hj = cos(z/j) e x + sin(z/ J ) e y , j — 1, 2, 3 , (A14) 

and equations (|A12| IA13|1 , after some simple algebra, assume the following form 

Ff = [sm 2 (uj - 0) - cos 2 (^- - 0) cos 2 9} , (A15) 
F* = [cos 9 sin2(z/ j - 0)] , j = 1,2,3 . (A16) 

Equations (JA15[ IA16J) show that the Fj~ beam-pattern functions are invariant under the 
transformation 9 — > 7r — 9, while the F„* change sign. These symmetry properties of the 
Fp x functions imply that the ZSS changes sign under the same transformation, an important 
property highlighted in Section II. 
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APPENDIX B: THE ZSS COEFFICIENTS a u a 2 , a 3 



In order to derive the expression of the ZSS in the Fourier domain we have introduced 
the generalized beam-pattern functions a + ,a x , j3 + ,j3 x , and 7+,7 x (see equations Hoi Hoi 
IT7j) . while the corresponding expression in the time-domain ( equation I27j) was written in 
terms of eighty-one amplitudes, Aj k \ and corresponding time-delays, , j = 1, 2, 3 , k = 
1, 2, 27. In this appendix we provide the analytic expressions of these functions. 

Using the identities /ij = k ■ pi and 

i 

Wt = k ■ hi = — (/i 2 - /X3) > ( B1 ) 

i 

w 2 = k ■ h 2 = — (/i 3 - Mi) > (B2) 

i 

w 3 = k -n 3 = — (/ii - H2) , (B3) 
the signal part of a(i) is written in the time domain as 

a 9W 



= [1 - w 2 ] [* 2 (t - 
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- [l+w 2 ][* 2 (t- 
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- V 3 (t-fiil 


')] 








L x 


-L 2 - 
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(B4) 

Since in the frequency domain the functions are equal to 

~ F+h + + F x h x 

*' = 2(1 - W ?) ' (B5) 
the response of the interferometric variable a, to a gravitational wave signal, becomes 

a(f) = a + (f,OA)h+ + aMe,<f>)hx • (B6) 

where the coefficients a +iX (f, 9, 0) may be written as 

a +tX (f,e,<f)) = axF x + ' x + a 2 F 2 + ' x + a 3 F+' x . (B7) 



In equation (|B7|) the functions ax, ct 2 , a 3 are given by the following expressions 

/ e i27r/(L 1 +L 2 +M2^) _ e i27r/(L 2 +M3<0) f {L 3 + nit) _ e t27r/(Li+£ 3 +M3*)) 

01 = 2|TT^) + 2(1^0 (B8) 
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( e i27r/(L 2 +M3<0 _ e i2nf( Ml e)\ f e i2T7f(L 1 +L 3 + f j, 3 £) _ e i2-Kf(L 1 +L 2 +L 3 + i i 1 t)\ 

"2 = ; ; + ; ; (B9) 

2(1 + W 2 ) 2(1 -w 2 ) K } 

/ e i27r/(Li+L 2 +i3+^i^) _ ^/(Lx+Z^+^h Ui2-Kf{^t) _ e i2nf(L :i +ti 2 e)\ 

a 3 = ; r + ; r (BIO) 

2(1 + w 3 ) 2(1 -w 3 ) 1 ; 

The numerical indices are again spacecraft indices, which can be permuted in order to obtain 
the corresponding functions for the other TDI combinations (3, and 7 (indices 1 — > 2 — > 3 — ► 
1, such that a — > (3 — ► 7 — > a). 

The 03 coefficient in the Zero Signal Solution is built from combinations of the pattern 
functions that determine the responses of a and j3 to the gravitational wave signal. Its 
analytic expression can now be rewritten as 

a 3 = a + (3 x - a x (3 + 

= (a 1 /3 2 -a 2 A)(F 1 + F 2 x -F 1 x F+) 
+ (a 3 ft - « 1 /3 3 )(F 3 + F 1 X - F 3 X F+) 

+ (a 2 (3 3 - aM{F+F* - F 2 X F 3 + ) . (Bll) 

After algebraically expanding the expression above we can rewrite 03 in the following 
form 

as = E 4 fc) e 2m/A 3 fe) , (B12) 
k=i 

where the are 27 different time shifts associated with the coefficient 03, and the 
are functions that only depend on the geometric angles (0,<f>). Since there are 27 unique 
time delays associated with each of the {a{\ coefficients, and because the gravitational wave 
signal enters in each TDI variable at 6 different times (the so called "six-pulse response"), 
it follows that in general the ZSS will display a total of "3 x (6 x 27) = 486 pulses". At the 
correct source location of course all these pulses add-up to zero! 

The analytic expressions for the 27 time delays entering into a 3 are given in the table 
below. The remaining 54 can be obtained from them by making a permutation over the 
spacecraft indices. 

Finally, the functions A 3 k \ introduced in equation (|B12|) . can be written in the following 
form 



^3 



Tg^itf - F X F 2 + ) + ^(F+F? - F*F+) + xiSi^ - F X F 3 + ) , (B13) 
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TABLE I: The twenty seven time delays entering into the function a 3 , which multiplies the Fourier 
transform of the measurement 7. The other fifty four time delays, entering into the functions a\ 
and a 2 , can be obtained by permuting, once and twice respectively, the spacecraft indices of the 
expressions below. 
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